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DIOCOTRON  INSTABILITY  OF  A  RELATIVISTIC  COAXIAL 
MULTI-RING  HOLLOW  ELECTRON  BEAM 

I.  Introduction 

There  has  been  considerable  interest  in  recent  years  in  the 
development  of  powerful  relativistic  electron  beams.1-3  Intense 
relativistic  hollow  electron  beams  have  been  used  in  the  laboratory 
recently  for  high-power  microwave  generation2,  in  electron  ring 
accelerators  and  autoaccelerators3.  .  .  etc. 

The  autoacceleration  process  is  a  collective  acceleration 
mechanism  for  generating  a  high  kinetic  energy  hollow  electron  beam. 
Such  beams  are  often  designed  to  have  approximately  uniform  electron 
density  within  an  annular  space  near  the  wall  of  a  cylindrical  drift 
tube,  with  no  electrons  outside  the  annulus.  Previous  theoretical 
studies  of  the  diocotron  instability  in  hollow  beams  have  assumed  a 
beam  profile  of  this  type. 

We  have  considered  a  somewhat  more  general  type  beam  profile, 
namely  a  hollow  coaxial  multi-ring  electron  beam.  Using  a  macroscopic 
cold  fluid  model,  the  diocotron  instability  which  characterizes  a  hol¬ 
low  multi-ring  electron  beam  has  been  investigated  for  a  broad  range 
of  beam  parameters  and  different  geometries.  Significantly  different 
results  exist  between  the  multi-ring  and  single-ring  hollow  beams. 
Description  of  these  differences  is  the  purpose  of  this  paper. 

The  basic  theory  and  assumptions  are  described  and  equilibrium 

properties  of  the  beam  are  examined  in  Section  II.  In  Section  III  we 

confine  our  study  to  systems  exhibiting  linear  behavior  from  which  an 

eigenvalue  equation  for  the  perturbed  potential  is  derived;  In  the 

case  of  a  square  radial  density  profile  a  closed  algebraic  dispersion 

relation  for  the  complex  eigenf requency  is  extracted.  The  dispersion 
Manuscript  submitted  June  4,  1980 


1 


« 


relation  is  solved  numerically  in  Section  IV  and  stability  properties 
are  investigated  in  detail.  Conclusions  are  drawn  in  Section  V. 

II.  Basic  theory 

We  start  by  considering  a  two-ring  electron  hollow  beam  with  a 
smooth  perfectly  conducting  wall  as  shown  in  Fig.  1.  Analysis  of 
dynamic  properties  is  based  on  a  macroscopic  cold  fluid  model  which  is 
idealized  in  that  the  flow  is  laminar  and  there  is  no  variation  in  the 
axial  direction.4  We  assume  a  system  infinite  in  the  axial  direction, 
and  a  strong  uniform  background  magnetic  field  which  prevents  the  beam 
from  spreading.  We  consider  a  cylindrically  symmetrical  electron  beam 
containing  n  electrons  per  unit  volume  moving  along  a  strong  axial 
magnetic  field  in  a  conducting  drift  tube  at  low  pressure  so  that  the 
beam  is  neither  current  nor  charge  neutralized.  Two  components  are 
assumed,  ions  with  charge  q  have  no  component  of  velocity  along  the 
beam  and  electrons  move  with  velocity  Pc  e z  which  is  assumed  large  com¬ 
pared  with  the  transverse  velocity.  We  allow  partial  neutralization 
by  a  fraction  f  of  charges  of  opposite  sign  trapped  in  the  beam.  The 
neutralization  fraction  f  denoting  the  ratio  of  neutralizing  charge  to 
beam  charge  is  assumed  uniform  across  the  beam.  By  virtue  of  these 
assumptions  the  electrons  are  described  in  cylindrical  geometry 
(r,  6,  z)  as  a  macroscopic  cold  fluid  immersed  in  a  uniform  axial 
magnetic  field  Bq  ez  with  both  radial  space-charge  and  azimuthal  self- 
magnetic  fields  included.  The  continuity  equation  and  the  equation  of 
motion  for  the  electron  fluid  can  be  expressed  in  the  relativistic 
form  as 
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Sn/St  +  V • (nV)  =  0 


(1) 


(d/St  +  V*V)  ymV  =  q(E  +  V  x  B)  (2) 

where  n(x,t)  and  V(x, t)  are  the  density  and  mean  velocity  and  E(jc,  t) 
and  JKjc,  t)  are  the  electric  and  magnetic  fields  respectively, 
q  and  m  are  the  charge  and  rest  mass  of  the  electron  and  y  =  (l-g2)-"5 
and  P  =  v^/c  are  the  standard  relativistic  quantities  and  c  is  the 
velocity  of  light  in  vacuum.  The  system  can  be  closed  by  adding 
Poisson's  equation  and  Ampere's  law  as  shown  below  respectively. 

V*E  =  eQ  ^l-f)  q  n  (3) 

V  x  B  =  p  qgcne  +  p  e  SE/St  (4) 

—  o  zoo  — 

where  p^  and  are  permeability  and  permittivity  of  free  space. 

The  equilibrium  state  (S/St  =  0)  is  azimuthally  symmetric 
(S/SO  *  0  and  S/Sz  =  0)  and  is  characterized  by  electron  density  n(r) 
and  azimuthal  electron  fluid  velocity  Vg  e^ .  The  deviation  from 
charge  neutrality  produces  a  radial  electric  field  that  influences  the 
azimuthal  motion  of  the  electron  fluid.  In  the  case  of  a  sharp¬ 
boundary  equilibrium  in  which  the  electrons  have  a  double  rectangular 

density  profile  as  shown  in  Fig.  2,  where  r  *  R  is  the  radial  location 

c 

of  a  grounded  conducting  wall,  the  self-generated  radial  electric  and 
azimuthal  fields  can  be  obtained  by  integrating  equations  (3)  and  (4). 
Thus, 
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Er(r)  =  (l-f)p_1B0(r) 


(r2  -  Rx2)/r 


=  ~  q  n(l-f) j  (R22  '  Rl2)/r 


(R22  "  Ri2  +  r2  -  R32)/r 

(R  2  -R2+R2-R  2 ) /r 
■^143 


Rl  <  r  <  R2 
R2  <  r  <  R3 
R3  <  r  <  \ 


R,  <  r  <  R 
4  c 


The  radial  electric  field  arising  from  the  space  charge  has  been 
reduced  by  a  factor  of  ( 1— f )  because  the  effect  of  partial  neutraliza¬ 
tion  by  ions.  It  follows  from  Eq.  (2)  that  equilibrium  force  balance 
in  the  radial  direction  can  be  expressed  as 


,(Er  +  VeBo  -  3c  B0> 


Eq.  (6)  is  simply  a  statement  of  radial  force  balance  of  centrifugal, 
magnetic  and  electric  forces  on  an  electron  fluid  element.  The  self- 
magnetic  field  produces  a  force  towards  the  axis  which  is  weaker  than 


the  outward  radial  electrostatic  force.  The  balance  among  electric, 
centrifugal  and  magnetic  forces  gives  the  angular  velocity  uu^Cr)  of  an 


electron  fluid  element  in  slow  rotational  equilibrium 


,  v  _  V6  ^pb  l-y2f 

Vr)  = 


(rz  -  R12)/r2  Rj  <  r  <  R2 

(R  2  -  R  2  +  r2  -  R  2)/r2  R,  r  <■  R, 


Rj  <  r  <  R2 


where  u)^  and  ui  ^  are  the  electron  cyclotron  and  plasma  frequencies 
respectively,  uu^  can  be  permitted  to  depend  on  r,  giving  sheared  cold 
fluid  rotation  rather  than  rigid  rotation.  Laminar  flow  and  the 
assumption  of  azimuthal  symmetry  together  imply  that  individual  charges 
moves  in  helices  of  constant  radius. 
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III.  Stability  analysis 

We  assume  all  the  perturbed  quantities  satisfy  the  conditions 
d[  ]/dz  =  0,  5[  ]/dt  =  iui[  ],  and  d[  ]/d0  =  — i [  ]  with  Im(uu)  <  0, 

where  l  is  the  azimuthal  harmonic  number.  After  Fourier  decomposing 
the  f luid-Maxwell  equations  (1)  to  (4),  it  is  straightforward  to  show 
that  in  the  strong  magnetic  field  regime  uucb  »  uu^  the  eigenvalue 
equation  for  the  perturbed  field  has  the  form 


(r  lr c  If  -  ?)6<t(t)  - 


-2i6cp(r)[6(r-R1)  -  6(r-R2)  +  6(r-R3)  -  6(r-R4)] 


r[  — 


tob(r) 


] 


(8) 


where  the  perturbed  potential  6cp(r)  =  6<|>  (r )  -  g  6  Az(r),  <f>  and  A  are 
the  scalar  and  vector  potentials  of  the  electromagnetic  field,  the  dio- 
cotron  frequency  is  defined  by  =  UJpb2/2y2UJcb,  and  ^b^  was  8iven  in 
(7).  The  right-hand  side  of  the  eigenvalue  Eq.  (8)  is  equal  to  zero 
except  at  the  surface  of  the  beam.  Moreover,  the  eigenfunction  6cp(r) 
satisfies  the  vacuum  Poisson  equation  except  at  r  =  Rj ,  ,  R3  and  R4  , 

therefore  the  piece-wise  solution  for  the  homogeneous  equation  can  be 
expressed  as 


6<p(r) 


l  .  ,  -i 

a  r  +  b  r 


l  ,  ,  -t 

c  r  +  d  r 

l  .  e  -1 
e  r  +  f  r 

,  -l 
g  *  +  h  r 

.  «...  -I 

i  r  +  j  r 


0  <  r  <  Rj 

R1  <  r  <  r2 
R2  <  r  <  r3 

r3  <  r  <  R4 

R  <  r  <  R 
4  C 


(9) 
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The  ten  coefficients  are  functions  of  Rp  RJt  R?  .  .  .  r 
to  be  determined  by  the  boundary  conditions,  the  requirement  implied  by 
delta  functions.  The  eigenfunction  is  continuous  at  each  boundary  and 
vanishes  both  at  r  =  0  and  r  =  R^.  The  effect  of  the  delta  function 
can  be  considered  by  multiplying  both  sides  by  r  and  integrating  over 
the  infinitesimal  interval  from  r(l-e)  to  r(l  +  e)  with  e  -*■  0  in  the 
vicinity  of  r  =  Rj,  R2»  and  R^  respectively.  Therefore  we  have  the 
following  dispersion  relation  given  by  10  equations  with  10  unknowns 


b  =  0 

a  Ei*  +  b  *r‘  -<=«,*  +  <* 

C  R  1  +  d  R  ~l  -  e  R  1  +  f  R  ~l 
2  2  2  2 

—  pi  r  n  D  1  ,  .  D  ~  2 

e  K2  t  1*2  ~  C  R^  +  d  R2 

eR3l  +  f  V*  “  8  R3£  +  h  R3~£ 

«  R3£  -  h  R3_£  ~  e  R3£  +  f  R3~£ 

8  \l  +  h  R4“*  -  i  R4£  +  j  R4~* 

i  R4A  -  j  R  ~l  -  g  R  1  +  h  R  ~l 

4  4  4  4 

iRJ  +  J  R  ~l  -  0 

C  C 


-2  (a  RjA  +  b 

UD 


2(c  R2£ 

a; _ 

%  ' 


+  d  R2_£) 


-2(e  R^  +  f  R3_i) 


2(g  R 4  +  h  R  "*) 


0) 

—  “  Y 
lbD  3 


(10) 
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where  Y  =  <l-y2f)  (R.,2  -  R12)/R22 

Y2  =  (l-y2f )  (R22  -  R^/R.,2 

Y3  =  (l-y2f)  (R42  -  R32  +  R ^  ~  Ri2>/R42 

The  determinant  of  the  10  linear  equations  gives  the  dispersion  relation 
of  the  form 


C1  =  R1C  -  R2c  +  R3c  -  V  “  Y1  -  Y2  -  Y3 

C2  "  R12  -  R13  +  R14  +  R23  "  R24  +  R34  +  Y1  "  Y2  +  Y3  '  2  +  V2 

+  V3  +  Y1Y3  -  Rlc(1  +  Y1  +  Y2  +  V  +  R2c(1  +  Y2  +  Y3> 


'  R3c(1  +  Y!  +  V  +  \c(1  +  Y1  +  Y2> 

C3  "  Y1  +  Y2  +  Y3  -  2Y1Y3  “  Y1Y2Y3"  R12<Y2  +  V  +  R13(Y1  +  V 

-  R14<Y1  +  Y2>  -  R23Y3  +  R2  4Y2  "  R34Y1  ‘  ~  2Y2  "  Y1Y2 

-  Y1Y3  -  Y2Y3  -  R43>  +  R2C<X  -  2Y2  “  Y2Y3  “  R43> 

-  R  (1  -  2Y  -  Y  Y  -  R,  . )  +  R  (1  -  2Y  -  Y  Y  -  R  ) 

3c  1  13  12'  4C  1  12  12 

c4  =  (Yj  -  1)(Y2  +  1)(Y3  -  1)  +  Rlc[(l  +  Yj)(l  -  Y2)(l  +  Y3) 

-  (i  +  yx)r43]  -  r2c[(i  -  y2)(i  +  y3)  -  r43] 


-  R3c(1  +  V<Y1  '  1  +  R12>  +  R4c(1  +  Y2)(Y1  -  1  +  R12> 

-  R12[(Y2  +  1)(Y3  -  1)  -  R34]  +  R13(l  +  Y3)(l  -  Y3) 

-  R14d  -  Y2H1  +  V  -  R23a  -  V  +  R24<!  -  V  +  R34<Y1  - 

2  0  0% 

where  R  ^  denoting  the  ratio  of  R^  to  R^  If  we  let  R2  »  R3  then 
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Eq.  (11)  gives  the  results  of  mono-ring  hollow  beam  case  as  shown  in 
Uhm  and  Siambis^.  By  the  methods  analogous  to  those  represented,  it  is 
straightforward  to  extend  the  case  to  multi -ring  beams. 

IV.  Numerical  results 

It  is  concluded  in  Section  III  that  the  dispersion  relation  has 
quadratic  form  for  the  single-ring  case  and  quartic  for  two  rings,  and 
so  on;  i.e.,  the  order  of  the  polynomial  equation  is  twice  the  number 
of  rings.  The  dispersion  relation  (11)  is  solved  numerically  for  the 
complex  eigenf requency  uu  s  uu^  +  iUL  with  real  oscillation  frequency  ui 
and  growth  rate  for  the  unstable  mode.  One  important  feature  of 
Eq.  (11)  is  that  the  complex  eigenf requency  is  linearly  proportional  to 
the  diocotron  frequency.  Consequently  the  applied  magnetic  field 
strongly  reduces  the  growth  rate  for  fixed  beam  density.  As  a  result, 
it  is  more  instructive  to  keep  the  beam  and  plasma  parameters  fixed  and 
study  the  growth  rates  for  different  geometries.  Fig.  3  shows  the 
growth  rate  for  a  two-ring  beam  versus  gap-length  d  of  the  rings  for 
different  mode  numbers  £  and  specified  values  of  f.  The  gap-length  d  is 
defined  as  the  distance  between  R2  and  while  we  fix  the  position  Rj 
and  R^  and  keep  the  beamwidth  of  each  ring  identical.  The  total  current 
of  the  beam  has  been  held  constant  for  various  geometries  so  that  the 
growth  rates  are  evaluated  at  the  same  beam  energy.  As  we  can  expect 
from  Eq.  (10)  the  growth  rate  has  a  strong  dependence  on  y2f.  For 
simplicity  we  specify  y  =  2  and  show  values  of  f  on  each  curve.  Note 
that  the  zero  gap  case  is  equivalent  to  a  single-ring  which  is  always 
stable  for  mode  £  =  1  as  observed  by  Uhm  and  Siambis&.  Therefore, 
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beams  with  various  gap-length  will  destabilize  the  i  =  1  mode  at 
least  for  .04  >  F  >  .005  as  can  be  seen  easily  from  Fig.  3.  Note  that 
there  is  no  instability  for  the  Si  =  1  mode  when  f  =  0.  The  growth  rate 
is  higher  for  larger  mode  numbers  ;  For  l  £  2  the  growth  rate  remains 
almost  constant  for  f  =  0  while  the  growth  rate  decreases  rapidly  as 
the  fractional  charge  neutralization  f  is  increased,  i.e.,  the 
unstable  modes  can  be  stabilized  easily  by  a  small  fractional  charge 
neutralization  with  appropriate  gap-length.  By  the  way,  the  real  fre¬ 
quency  for  the  unstable  modes  shown  in  Fig.  3  is  larger  for  higher  l 
but  does  not  exceed  .6  uip.  It  has  been  demonstrated  that  the  growth 
rate  of  the  instability  exhibits  a  sensitive  dependence  on  f. 
Neutralization  is  commonly  produced  from  the  residual  gas  in  the 
apparatus.  It  is  difficult  to  give  a  universal  theory8  for  f  which 
depends  very  much  on  such  factors  as  the  energy  of  the  beam  and  the 
composition  and  pressure  of  the  residual  gases,  ionization  cross 
section,  the  energy  of  the  secondary  electrons  and  whether  the  ion  can 
escape  from  the  end  of  the  beam. 

Another  interesting  feature  is  that  for  higher  modes  the  double 
region  in  gap-length  for  instability  disappears  from  Fig.  3  to  Fig.  4. 
Generally,  the  unstable  mode  for  a  single  ring  can  be  stabilized  by 
increasing  the  gap-length  of  a  two-ring  beam  in  this  kind  of  geometry. 
Next,  we  want  to  demonstrate  the  effect  of  the  total  beamwidth 
(R4  -  Rj)  on  the  diocotron  instability  which  is  illustrated  in  Fig.  5 
and  Fig.  6.  The  l  =  1  mode  is  unstable  as  before  when  the  single-ring 
beam  becomes  a  two-ring  beam  as  shown  in  Fig.  5.  However,  the  funda¬ 
mental  mode  can  no  longer  occur  if  the  total  beamwidth  has  expanded  to 
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the  one  as  shown  in  Fig.  6.  Clearly,  the  i  =  1  mode  can  be 
destabilized  by  moving  two  rings  together.  For  l  i  2,  the  growth  rate 
increases  with  respect  to  gap-length  when  f  =  0;  Nevertheless,  the 
growth  rate  for  small  f  increases  with  respect  to  gap-length  first  and 
then  decreases  sharply. 

Finally,  it  is  straightforward  to  extend  the  analysis  to  a  three- 
ring  geometry  as  shown  in  Fig.  7.  For  easy  comparison  between  the  two- 
ring  and  three-ring  cases,  we  plot  the  growth  rate  versus  gap-length 
d  =  (R^  -  R^)  while  fixing  the  positions  R^  and  R6  and  keeping  the  beam- 
width  of  each  ring  identical,  the  total  current  is  also  fixed.  The 
effect  of  inserting  the  center  ring  can  be  seen  easily  while  comparing 
Fig.  8  to  Fig.  3.  The  unstable  modes  in  two-ring  geometry  have 
been  retained  in  the  three-ring  case,  but  also  the  growth  rate  has  been 
enhanced  by  inserting  the  third  ring  in  the  middle,  which  may  be  caused 
by  the  strong  coupling  among  the  self-fields  of  the  beams.  For  higher 
mode  numbers,  unstable  modes  can  be  found  in  three  separate  regions  in 
gap-length  domain.  If  we  move  the  outer  edge  of  the  beam  away  from  the 
wall.  Fig.  9  gives  the  geometry  which  can  be  compared  to  Fig.  6.  It 
comes  as  no  surprise  that  the  l  =  1  mode  becomes  unstable  again.  For 
higher  modes,  the  triple  region  for  instability  is  not  obviously  seen; 
The  conducting  wall  plays  an  important  role  in  the  diocotron 
instability. 

V.  Conclusion 

We  have  formulated  a  f luid-Maxwell  theory  of  the  diocotron 
instability  in  an  infinitely  long  relativistic  electron  beam  propagating 
parallel  to  a  uniform  applied  axial  magnetic  field.  In  beams  with 
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self-fields  it  frequently  permits  simple  models  which  illustrate  many 
of  the  essential  features  of  more  realistic  types  of  beam.  The  growth 
rate  has  been  calculated  with  special  emphasis  put  on  displaying 
results  as  a  function  of  conducting  wall  geometry,  i.e.,  conducting  wall 
location  relative  to  the  beams  position.  The  results  show  the  strong 
influence  of  neutralization  fraction  f,  relative  position  of  the  rings 
with  conducting  wall  and  gap-length  of  the  rings  on  the  diocotron 
instability.  It  seems  the  instability  has  been  enhanced  by  the  multi¬ 
ring  geometry  in  a  rather  complicated  manner.  In  short,  the  l  =  1  mode 
which  is  always  stable  in  the  single-ring  case  can  become  unstable  in 
the  multi-ring  geometries.  However,  in  the  two-ring  geometry,  the 
Si  =  1  mode  diocotron  instability  can  be  avoided  by  either  moving  two 
hollow  beams  away  from  the  wall  or  spreading  the  two  rings  farther 
apart.  For  the  H  £  2  modes,  the  growth  rate  can  be  either  enhanced 
(e.g.  Fig.  5  and  Fig.  6)  or  reduced  (e.g.  Fig.  3  and  Fig.  4)  depending 
on  the  position  and  gap-length  of  the  rings.  In  the  three-ring  case, 
all  the  unstable  modes  occurring  in  the  twO'-ring  geometries  have  heen 
retained  and  the  growth  rates  have  been  enhanced  at  the  same  gap-length. 
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